ESSSA Group 


Cx \ 

NASA! 


A Bright Lunar Impact Flash 
Linked to the Virginid Meteor Complex 

m- 


D. E. Moser 

Jacobs, ESSSA Group, Meteoroid Environment Office, NASA Marshall Space Flight Center 


R. M. Suggs, R. J. Suggs 

NASA, Meteoroid Environment Office, NASA Marshall Space Flight Center 


Stanford Meteor Environment and Effects Workshop 
Stanford University, California, 14-16 July 2015 



On 17 March 2013 at 03:50:54 UTC, NASA detected a bright impact 
flash on the Moon caused by a meteoroid impacting the lunar surface. 


There was meteor activity in Earth’s atmosphere the same night from 
the Virginid Meteor Complex. 


The impact crater associated with the impact flash was found and 
imaged by Lunar Reconnaissance Orbiter (LRO). 
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9-year observing program 



Goal: Monitor the Moon for impact flashes produced by meteoroids striking the lunar surface. 
Determine meteoroid flux in 10’s g - kg size range. 



Observation from MSFC Field of View 

• Two 0.35m telescopes simultaneously • FOV covers approx. 20 arcmin 

• Black & white CCD video cameras • 4* 10 6 km 2 on the leading or trailing edge 

• Interleaved 30fps video digitized, recorded • Observing when illumination 10-50% 

• Video analyzed with custom software • Maximum 10 observing nights/month 
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330+ Lunar impacts observed 
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For more info: 

Suggs et al. (2014) 
Suggs et al. (2011) 
Suggs et al. (2007) 
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Typical impact flashes 


15 Dec 2006 
09:17:39.336 
33 ms 


m R = 7.4 


□ 


0.09 kg 

Geminid (35 km/s) 


3 Nov 2008 
00:11:06.144 
100 ms 


m R = 7.7 


□ 


0.1 kg 

S. Taurid (27 km/s) 






17 Nov 2006 
10:56:34.820 
66 ms 

m R = 7.0 Q 

0.03 kg 

Leonid (7 1 km/s) 


22 Apr 2007 
03:12:24.372 
133 ms 


m R = 6.7 


0.08 kg 

Lyrid (49 km/s) 
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iSSfSIS Atypical flash on 17 March 2013 


17 Mar 2013 
03:50:54.312 
1.03 s 

m R = 3.0 ± 0.4 
16^° kg 
Virginid 
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Flash info 

Detected with two 
0.35 m telescopes 

Watec 902H2 Ultimate 
monochrome CCD 
cameras 

— Manual gain control 
a - gamma = 0.45 

CO 

cd 

| Saturated — ► needed 

O) 

1 saturation correction! 

■O 

DO 

= - Photometry performed 

1 using comparison stars 

| (Suggs et al. 2014) 

DjO 

c 

^ - 2D elliptical Gaussian 

I fit to unsaturated wings 

S ofPSF 

-Q 

o 

Peak m R = 3.0 ± 0.4 
Lum energy = 7.1 t^x 10 6 J 


Jacobs, ESSSA Group/MEO/D.E. Moser 


SMEE Workshop, 14-16 July 2015 


6 



JACOBS 

ESSSA Group 


02:17:00 
33 ms 



03:14:21 
66 ms 



Increased lunar activity 
on 17 Mar 2013 



o3 

I £r 


02:33:18 
33 ms 


03:50:54 
1.03 s 
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Observed by A. Kingery & R.M. Suggs; detected by R.J. Suggs Obsen/ecnDy A. Kingery & R.M. Suggs; detected by R.J. Suggs 




iSSfSS Meteor activity on 17 Mar 2013 
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19 meteors were observed by NASA & SOMN 
all-sky meteor cameras on 17 Mar 2013 
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Meteor shower activity 
on 17 Mar 2013 (SOMN & NASA) 



19 meteors seen on 17 Mar 2013 Cluster of 5 seen on 17 Mar 2013 



RA (J2000.0, deg) 


Geocentric meteor radiants color-coded by 
speed with a tight cluster of 5 with: 


\ 

x $ 

Vv • j 


Orbits of the cluster of 5 were very similar 
with the following average orbital elements: 

meteoroids 

a (AU) 2.25 ±0.17 
e 0.79 ± 0.02 

i(°) 5.26 ±1.02 

to (°) 280.32 ±2.11 

Q (°) 356.65 ± 0.07 

q (AU) 0.48 ± 0.02 
Q (AU) 4.0 ±0.3 


meteors 

a g (°) 184.1 ±1.0 

8 g (°) 4.4 ± 0.9 

v g (km/s) 25.6 ± 0.8 

^sun (°) 356.6 


Tj 3.1 ±0.2 — > Indicates -asteroidal body 

SMEE Workshop, 14-16 July 2015 
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Jacobs Meteor shower activity 
EsssAGroup on 17 Mar 2013 (SOMN & NASA) 


19 meteors seen on 17 Mar 2013 Cluster of 5 seen on 17 Mar 2013 



Geocentric meteor radiants color-coded by Orbits of the cluster of 5 were very similar 

speed with a tight cluster of 5 with: with the following average orbital elements: 


Virginid Complex 


meteoroids 

NVI 

EVI 

(rad at k=356.6) 

a (AU) 

2.25 ±0.17 

1.69 

2.82 

meteors NVI 1 EVI 2 

e 

0.79 ± 0.02 

0.71 

0.86 

a g (°) 184.1 ±1.0 183.1 181.0 

i(°) 

5.26 ± 1.02 

3.7 

5.2 

8 g (°) 4.4 ±0.9 2.3 4.7 

<0(°) 

280.32 ±2.11 

282.4 

285.8 

v g (km/s) 25.6 ±0.8 23.0 28.9 

Q(°) 

356.65 ± 0.07 

358.0 

355.1 

X sun (°) 356.6 354 354 

q (AU) 

0.48 ± 0.02 

0.496 

0.40 


Q(AU) 

4.0 ±0.3 

2.89 

5.25 

'Sekanina (1973), 2 Whipple (1957) 

Tj 

3.1 ±0.2 

Indicates -asteroidal body 
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Favorable Virginid 
radiant geometry 



iflMmpactFlash 
. subradiant 


Pink indicates the portion of the moon visible to the radiant. 
Impact angle -56° from horizontal. 
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Putting it all together 



All- sky meteor cameras detected a deeply 
penetrating cluster of 5 meteors on 17 March. 

Radiant and orbital elements consistent with the 
Virginid Meteor Complex (EVI/NVI). 

Impact flash rate increased to 1 every 0.87 hours 
on 17 March. (4 impacts in 3.5 hours) 


Impact Constraints 

Assume impact flash was part of 
Virginid Meteor Complex 
v g = 25.6 km/s 
% = 56 ° 

^ Asteroidal? (Tj = 3.1 ± 0.2) 
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Mapping the impact location Sm 



ArcMap was used to selenoreference the 
lunar impact following a newly 
developed selenolocation workflow. 


Using the intensity-weighted center of the flash 


Euler 

MARE IMBRIUM 


+ 

Pytheas 

.. ^ ^ 

Draper 

i . 1 i i i"i i 

•ft | 1 


Final predicted crater position 
20°. 684 ± 0.258 N, 24°. 228 ± 0.288 W 
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Impact crater found by LRO! 

Robinson et al. (2014) 



March 1 ? h Impact 

|»' IS ■> » A 





Image from Robinson (2013) NASA/GSFC/Arizona State University 


50m 


Features 

• Fresh, bright ejecta 

• Circular crater 

• Asymmetrical ray pattern 


Crater info 

• Rim-to-rim diameter = 18 m 

• Inner diameter = 15 m 

• Depth ~ 5 m 

Actual crater location 
•20.7135° N, 24.3302° W 

Impact Constraints 
^ Circular crater, impact 
angle constrained 0 h >15° 
^ Ejecta gives no azimuth 

constraint (Robinson, personal comm. 2014a) 
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Comparison of geolocation 
results to obs crater location 
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Transient crater diameter estimates 

Assumptions: Virginid v gfoc =25.7 km/s, 0 h = 56°; p t = 1500 kg/m 3 (regolith) 



Model 

Lum eff. r| 

KE 

xl0 9 (J) 

Mass 

(kg) 

Pp 

(kg/m 3 ) 

Dcaic 

(m) 

Dobs 

(m) 

% Err 

Gault's crater 
scaling law 

(Gault 1974) 

5 X 10' 4 

(Bouley et al. 2012) 

14 

[9.4,22] 

42 

[28,66] 

1800 

18.5 [16.5,21.1] 

15 

23% 

3000 

20.2 [18.0,23.0] 

15 

35% 

1.3X10- 3 

(Moser etal. 2011) 

5.4 

[3.6, 8.4] 

16 

[11,26] 

1800 

14.1 [12.5,16.0] 

15 

6% 

3000 

15.3 [13.6,17.4] 

15 

2% 

Holsapple's 

online 

calculator 

(Holsapple 1993) 

5 X 10" 4 

(Bouley et al. 2012) 

14 

[9.4,22] 

42 

[28,66] 

1800 

12.2 [10.9,13.8] 

15 

19% 

3000 

12.5 [11.1,14.2] 

15 

17% 

1.3X10 3 

(Moser et al. 2011) 

5.4 

[3.6, 8.4] 

16 

[11,26] 

1800 

9.3 [8.3,10.5] 

15 

38% 

3000 

9.5 [8.5,10.8] 

15 

37% 


Two example values of r\ from the literature yield large ranges for KE and mass. 
Consequently, model results are highly dependent on luminous efficiency r|. 


Assuming a velocity dependent r\ = 1.3xl0 3 , these model 
results are consistent with the observed crater diameters. 

D ca i c = 8-18 m transient crater D obs = 15 m inner (‘transient’) 

D calc = 10-23 m rim-to-rim D obs = 18 m rim-to-rim 
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Impact summary 
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Date of impact: 
Duration of impact: 
Corrected flash peak R magnitude: 
Luminous energy generated by impact: 
Estimated kinetic energy of impactor: 
Estimated mass of impactor: 
Estimated diameter of impactor: 

Crater diameter: 
Crater location: 
Possible meteor shower association: 


17 March 2013 3:50:54 UTC 
1.03 s 

3.0 ±0.4 
7 . 1 S 4 X 10 6 J 

5.4-u x 10 9 J = 1.3 tons of TNT (assuming r| = 1.: 

16? kg (assuming v = 25.7 km/s) 

22 ± 3 cm (assuming p p = 3000 kg/m 3 ) 

18 m rim-to-rim, 15 m inner (‘transient’) 
20.7135° N, 24.3302° W 

Virginid Meteor Complex 


X 10-3 ) 
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Backup Slides 



Luminous efficiency: LE a = HaKE 
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D t = 15 m from crater measurements "1 Gault’s crater scaling law (Gault 1974) rearranges to give 

LE A = 7. 1 x 10 6 J from flash measurements J jj a vs 0 h without assuming impact speed. 

nx = LE a /KE = LE A / (4.0 APp 0167 A 0 ' 5 sin- 1/3 A) 1/0 - 29 



Projectile densities 

- — * pp= 1000 kg/m J 

— p p = l S00 kg/m-'' 

*— pp = 2B00 kg/m 3 

pP = 3000 kg/m-"' 

— pP = 330(1 kg/m 3 


Typical values of r/ A derived from lunar regolith range from 2 ± 1 x 1 O' 4 to 2 ± 1 x 10' 3 . 
Assuming association with the Virginids, 6 h = 56° and 7.52.5 x 10' 4 < iJ A < 1.52!* x 10' 3 . 
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Crater scaling laws 
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Gault's scaling law (Gault 1974) for D < 100 m 
D = 0.25pp 0167 p t _0 - 5 £ 0 - 29 (sin6i ft ) 1/3 

D = transient crater diameter 

p p = projectile density 

p t = target density 

E = kinetic energy of projectile 

6 h = impact angle measured wrt horizontal 


Holsapple crater calculator 

(http://keith.aa.washington.edu/craterdata/scaling/index.htm) 


Crater Sizes from Explosions or Impacts 



1 1 r'JL* i Mfr'fnvi | 

Eftcte unt 'hie mu men ymstm 

t» or N "♦ 
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• Transient crater diameter • Crater rim diameter 

• Measured at the height of the • Measured from rim to rim 

pre-existing surface • Aka outer diameter 

• Aka inner diameter, apparent 
diameter 

Holsapple assumes D rim = 1.3 D similar to Melosh (1989) 
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Ejecta distribution 

after Robinson et al. (2014) 




Ejecta in multiple reflectance zones 

High reflectance zone 10-20 m SW, <10 m NE 
Low reflectance zone 50 m WSW, 80 m ENE 

-300 m rough semicircle 
Low reflectance zone -1 km centered in NE 
248 circ/irreg splotches within 30 km See Robinso1 

O r for mor 


See Robinson et al. (2014) 
for more details 


Impact Constraints 

^ Circular crater, impact angle constrained >15° 

HRZ - impact possible from SE or NW 
^ LRZ - impact possible from SW 

.*. no azimuth constraint (Robinson, personal comm. 2014a) 


An impact from the SW is consistent with an 
impactor from the Virginid Meteor Complex. 
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Peak Impact Flash Magnitudes 



40 1 1 1 1 1 1 r 



ft Magnitude 

Mar 17 impact 
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All sky heights & masses 




Number of meteors 

Cluster of 5 penetrated fairly deeply into the atmosphere. 
Moderately large, 0.003 kg to 0. 1 kg masses 
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Crater associated with 
Spanish flash on 11 Sep 2013 




Image from Robinson (2014b) 


NASA/GSFC/Arizona State University 


Crater info 

• Rim-to-rim diameter = 34 m 

• Ejecta effects extend more 
than 500 m from the crater in 
all directions 
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